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Abstract 

We study the effect of thermal fluctuations in a recently proposed protocol 
for transmission of unknown quantum states through quantum spin chains. We 
develop a low temperature expansion for general spin chains. We then apply this 
formalism to study exactly thermal effects on short spin chains of four spins. We 
show that optimal times for extraction of output states are almost independent 
of the temperature which lowers only the fidelity of the channel. Moreover we 
show that thermal effects are smaller in the anti-ferromagnetic chains than the 
ferromagnetic ones. 
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1 Introduction 



One of the basic ingredients of quantum communication is the transport of a known 
or unknown quantum state from one point to another P^. Recently it has been shown 
that quantum spin chains can act as channels for the efficient [H El 01 or perfect |B] 
transfer of quantum states. Of particular interest to us here is the works reported 
in 1210] • It has been shown that this scheme may be used for linking several small 
quantum processers in large scale quantum computing. In this scheme, an A^-site 
ferromagnetic Heisenberg spin chain placed in a magnetic field, plays the role of the 
quantum channel. The spins of the channel are numbered from 1 to A^. It is assumed 
that this chain is in its ground state, |— , — , • • • — ), where |— ) denotes the state of a 
spin in the negative z direction and the magnetic field points to the positive z direc- 
tion. One then adds a spin or qubit to the left hand side of this channel labeled by 
zero, which is in an unknown state \^) (figure 1-a). The state of this qubit is to be 
transmitted with a high fidelity to the right hand side, by the natural time evolution 
of the chain. 

In this way one may circumvent a problem which exists in quantum computer imple- 
mentation namely the difficulty of switching on and off between spins |71|H]. 

The Hamiltonians governing the interaction of the spins in the channel and the 
full chain are respectively 

7V-1 TV 

and 



N-\ N 

H=-J^ai-ai+i+B^a,,u (2) 

where the subscript "c" on H stands for the channel. If at time t = the qubit 
is placed to the left of the chain then the evolution of the Heisenberg chain carries 
the state of this qubit to the rightmost spin where one extracts the state with a 
rather high fidelity, provided that one extracts the state at an optimal time. 
It has been shown in ^2] that one can transmit quantum states with a high fidelity 
ranging from F = 1 for = 4 to a value exceeding 0.9 for N = 7, 10, 11, 13 and 
14. The fidelity generally decreases with increasing the length of the channel exceeds 
the classical value of -F = 2/3 for A^ ~ 80 which holds for classical transmission of 
quantum states jH]. 

One can also use this channel for transmission of entanglement in the following way 
PI . One places two maximally entangled spins labeled 0' and to the left of the chain 
(figure 1-b) and evolution of the chain after a suitable lapse of time entangles the spin 
0' to the rightmost spin A'^ here it is assumed that only the spin is coupled to the 
channel. In this way two distant spins can be entangled which can later be used for 
implementation of other quantum protocols like teleportation jlUj . 
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Figure 1: a- An unknown quantum state ](/)) is placed at site of the chain and 
is transported to site N by the dynamics of the spin chain, b- The entanglement 
of a Bell state = "^(|0)0) + placed at sites 0' and develops into the 

entanglement of sites 0' and N. 2 



The formalism of P] requires that the spins of the channel be ideally aligned in 
the direction opposite to that of the magnetic field. This ideal situation is however 
achievable only at zero temperature or in very strong magnetic fields where thermal 
fluctuations are not large enough to populate excited states, i.e. (when B/kT >> 1). 
On the other hand increasing the magnetic field may lower the quality of the channel 
since a high magnetic field tends to align the spins and will generally dominate the 
interaction between the spins which is essential for the working of the channel. 
Moreover when one uses this channel once and extracts the state at the right hand 
site, the initial state of the channel turns into a mixed state. Before using the channel 
for another round of transmission the initial state of the channel should be restored, 
for example by cooling to low enough temperatures. It is plausible to assume that 
multiple uses of the channel may heat it up to temperatures in which not only the 
ground state but also some of the excited states are also populated. 

In view of these considerations it is desirable to study the effect of thermal fluc- 
tuations on such a quantum channel. This is the problem that we want to address in 
this paper. Thus we want to generalize the protocol of |5] to the case where the initial 
state of the channel is not the ground state but a thermal state given by a thermal 
density matrix. 

This enables us to to see the effect of the ambient temperature on the feasibility of 
the protocol, and the quantitative effect that temperature has on the fidelity of trans- 
mission of states and distribution of entanglement. 

We will derive a low temperature expansion through which we can study the effect 
of temperature to any desired degree of accuracy by keeping appropriate number of 
terms in the expansion. 

In a sequel to this paper we will study long chains of arbitrary number of spins at 
low temperatures. This study can be done only numerically. In this paper however 
we will study exactly thermal effects on a short chain of four spins. The advantage of 
studying this short chain is that we can obtain the spectrum completely and hence 
can compare the two cases of ferromagnetic and antiferromagnetic chains. 

The structure of this paper is as follows: In section [2 we set up the general 
formalism and will develop low temperature expansions for the expressions of fidelity 
and entanglement. In section 13] we derive general expressions for entanglement of 
endpoints of the chain at arbitrary temperatures. In section 0] we study exactly the 
specific example of a short chain of only 4 spins where we present our basic results in 
figures (HI to EI). 

2 Low temperature expansion of the fidelity 

We consider the interaction between the spins to be nearest neighbor and of Heisen- 
berg type. The spins also interact with an external magnetic field. We should em- 
phasize that much of what we derive in this section do not depend on a specific form 
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of the Hamiltonian. We assume that the initial state of the channel is 



Pth = —^ = 2_^ ^ (3) 

a 

where |a)'s and EaS are respectively the eigenstates and eigenenergies of the channel 
Hamiltonian (^. Here Z is the partition function of the channel Z := tr{e~^^''). 
At time t = we place the 0-th spin which is in an unknown pure state pQ = |</')(</'| to 
the left of this channel. The initial state of the whole chain will be given by po ® pth 
which evolves to 

p(t) = e-*^*(po®ft/^)e^^*, (4) 

where H is now the Hamiltonian of the full chain 0. 
The state of the A^-th spin after time t will be given by 

PN{t) = trj^{e-'"\po Pth)e''''), (5) 

where N means that we take the trace over all sites except the A^-th site. We can 
now derive an operator sum representation for the transformation of the state of the 
leftmost qubit p(0) := \(p){(p\ to the state of the rightmost qubit pNit)- To do this 
we use © to compute an element of pAr(t) as follows, where we use the index I to 
run over a complete set of states {|/)} for the Hilbert space of part of the total chain 
from site to site A — 1, and the index a is to run over the eigenstates of the channel 
Hamiltonian: 

{k\pN{t)\l) = Y.{I,k\e-'^\po®pth)e'''V,l) 
I 

= ^le"'''*!^' «)0' "K/'o » Pth)\m, f3){m, Pie'""']!, I) 

p — fiEa 

= Yl {lMe-'"%c^)p,,,m{m.aW"'\I,l)^^. (6) 

I,j,m,a 

Defining a collection of two by two matrices Mi^a with elements 

{k\Mi^^\j) := ^J^{I,k\e-'H%a) (7) 
vZ 



we find the following Kraus decomposition I12j for pjy (t) , 

PN{t) = YMic.PoM\^. (8) 

Note that in general the number of elements in this decomposition, i.e. the num- 
ber of matrices Mj ^ is huge. In fact this number is equal to the number of different 
choices for the pair of indices (/, a), which equals 2^^. In principle there are many 
operator sum representations for a superoperator and the number of kraus operators 
for a qubit can be reduced to 4 jlll I12j . however the present form has the advantage 
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that it is suitable for a low temperature expansion of the fidelity of the channel. Note 
also that even in the present form symmetry arguments highly restricts the number 
of nonzero matrices as we will see in the sequel. 



The fidelity of this out-state PN{t) with the in-state po is given by 

F = Y^ tr{poMiaPoMjj. (9) 

I,a 

We are interested in the fidelity averaged over all the initial input states, that is 

where the integral is taken over the surface of the Bloch sphere. This integral can 
further be simplified by using the following easily verified identity 

tr{poApoB) = tr{{A B)S{po po)) (11) 

where S is the swap operator with elements Sij^ki = ^u^jk- We now write po as 
Pq = ^(1 + n ■ a) and use the identity f nirijdft = ^6ij to arrive at the following 
identity 

I Po^Podn = ^{l + S), (12) 

which we use to rewrite F as 

F=^Y1 ® + S)). (13) 

I, a 

Using the facts that tr{{A ® B)S) = tr{AB) and Y.i,a ^j^Mja = 1 we find the 
final form of F as 

I,a 

Equation (|14|) is already in the form of a low temperature expansion for the average 
fidelity. The leading contribution comes from the ground state which we label as 
ao, the next to leading contribution comes from the first excited states and so on. 
Thus despite the huge number of matrices M/^q,, at low temperatures one can obtain a 
reasonably good value of the fidelity by using only the first few terms in the expansion. 



2.1 The zero temperature limit of the ferromagnetic chain 

In this limit only the ground state contributes to the expansion (|14|) . For the ferro- 
magnetic chain the ground state is 

|ao) = (15) 
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where all the spins are down. Thus we have 

F=1 + IY.\^''^^-o\'- (16) 

In H16|) it appears that a set of 2^ matrices Mj^ao contribute to the sum. However 
we show that by symmetry considerations one can reduce this number to only two. 
To this end we note that a general element of Mj^a^ can be written as follows: 

Mj,,„(i,j) = (/,i|e-^^*|j,-, ,-). (17) 

In view of the symmetry [H, Jz] = where Jz is the total spin in the z direction, the 
only nonzero matrices are those in which the index / is either (—,—,—,••• ,— ) for 
which we denote the corresponding matrix by Mq, or those in which only one spin 
is up (e.g. in the i-th position (— , —,•••—,+,—,••• , — )) for which we denote the 
corresponding matrix by Mj. 

Moreover the above mentioned symmetry requires that matrices be of the following 
form: 



Mi 







m+ 
m_ 



where 



and 



M. = ( ) ^ = 1..■N, (18) 



m+: = (-,-, ,+|e ^■^*| + ,-,-,--- ,-), 

m_ : = (-,-,••• -, -|e-*^*|-, (19) 



mi :={-,■ ■■-,+, |e-'^*| + , ")• (20) 



A simple calculation now shows that the following identity holds for the above 
types of matrices, regardless of their explicit form of matrix elements, 



N 



MipMj = MpM\ (21) 



i=l 

where 



M=( /^N , . )• (22) 



Thus the operator sum representation reduces to a sum with only two elements, i.e. 

PN (t) = M(,pqMI + MpoMt . (23) 
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Using (|16() and noting that the matrix M is traceless, we find the average fidehty 
at zero temperature: 

11 

F = - + -\mj^ + m^^ (24) 
6 D 

For the ferromagnetic chain, the state |— , — , • • • , — ) is also the ground state of the 
fuh Hamiltonian and thus by shifting the zero energy of the Hamiltonian to we can 
set m_ = 1. Thus we find 

F = ^ + i|l + m+|2 (25) 
in accordance with the result of ^2_. Note that our is denoted by fo,Nit) in [2]. 



3 Transfer of entanglement 

We now consider how temperature affects the distribution of a maximally entangled 
pair through the channel. Following |;2; we consider a maximally entangled pair of 
qubits in the state |$+) := "^(|0;1) + figure (1-b) this pair of qubits are 

labeled by 0' and 0. 

The evolution of the chain may transform the entanglement between the pair (0', 0) 
to the pair (0', A^), thus enabling us to transport entanglement between ions or any 
other realization of qubits over long distances. Note that the Hamiltonian only acts 
on the part of the chain form to N. It is assumed that the qubit 0' does not interact 
with the rest of the chain. After a time t, the density matrix of the pair (0',A^) is 
easily obtained thanks to the operator sum representation Q. We find 



P0',N{t) = ^(1 «) M/„)(|$+)($+)(l M 



t ) 

la) 



(26) 



I,a 



Using the fact that 



l</'+)('A4 



V 



1 1 
1 1 



(27) 



0/ 



and the explicit form of the matrix elements of Mj^a as given in Q we find after some 
manipulations the following low temperature expansion: 



PO'N 



z 



-P{a){t), 



(28) 



where each pa pertains to a level a of the spectrum of the channel and is given by 



w+it) z^it) 



\ 



(29) 



7 



where 



Za{t) 



i(l,a|(l+a^(t))|l,a) 
^(0,a|(l-a^(t))|0,a) 
i(l,a|(l-c7^(t))|l,a) 
l(0,a|(l + a^(t))|0,a) 
^(0,a|cT+(t)|l,a). 



(30) 



Here the operators cr^(t) are operators in the Heisenberg picture, i.e. (cr^(t) = 



Thus the concurrence will have the same dependence on the correlation function as 
in the thermal equilibrium state |14| I15j , only now the correlation functions are time 
dependent. 



4 Exact solution for a short chain 

An exact study of the thermal effects on a long chain with arbitrary number of spins 
is highly involved, since it requires a knowledge of all the energy eigenstates. On can 
study these chains only at low temperatures in which case only the ground state and 
the first excited states of the chain are populated. We will do this in a sequel to 
this paper. Here we study exactly a short spin chain of = 3. An exact study of 
a short chain has the advantage that one can compare the characteristically different 
behaviors of ferromagnetic and anti-ferromagnetic chains. We remind the reader that 
at zero temperature the channel relaxes to its ground state which for the ferromagnetic 
chain is a disentangled state of spins aligned with the magnetic field. This is the only 
case which has been studied in 0. 

So we consider a three site channel with hamiltonian given by 



The spectrum of this hamiltonian is easily obtained and it exist in the appendix. 
The Hamiltonian of the full chain is now 



Determination of the spectrum of this hamiltonian is facilitated by using the following 
symmetries, namely 




He = -J(S1 • S2 + S2 • S3) + B{si^ + S2z + S3z). 



(31) 



H = - J(so • Si + Si • S2 + S2 • S3) + B{soz + Su + S2z + ^Sz)- 



(32) 



[H, Jz] = [H, A] = [J^, A] = and a^^H{J, B) = H{J, -B)a^ 



(33) 



where A is the inversion operator 




(34) 
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Figure 2: (Color Online) The fidelity (F) between the output and the input states 
averaged over all input states, as a function of kT and time t, in a fixed magnetic 
field. In all the figures the fidelity F and the concurrence C are dimensionless and 
we are working in units where KT and t have no dimensions. In these units we have 
taken B = 1 and J = +1 for ferromagnetic and J = — 1 for anti- ferromagnetic chains. 

The spectrums of the total and the channel hamiltonian are derived in the appendix. 
By plugging these eigenstates and eigenvalues in equations Q and ((TH) one can 
determine the average fidelity of the output state with the input state. The fidelity 
is a complicated function of time, in fact it is a superposition of periodic functions 
with periods u>ij := |^ | . In 2 one extracts the output state only at certain 
times where the fidelity reaches a maximum. Since we want to focus on the effect of 
temperature, we fix the magnetic field and determine the average fidelity as a function 
of time and temperature. The results are shown in figure © and © for ferromagnetic 
and anti-ferromagnetic chains respectively. 

These curves show several interesting features. The first one is that the optimal 
time of extraction is almost independent of temperature, thus at any temperature one 
can tune the optimal time of extraction to be the same as that of the zero temperature. 
The only effect of temperature is that it decreases the fidelity. It is also seen that the 
optimal times when the fidelity reaches local maxima are the same for both types of 
chains. Moreover thermal fiuctuations have much less destructive effect on the fidelity 
in the anti-ferromagnetic chain as compared with the ferromagnetic chain. 
The difference between these two types of chains is more pronounced when we use 
them to transfer entanglement. We have used these two channels to transfer the 
maximal entanglement between the two spins 0' and at the left hand side of the 
chain Q to entanglement of the end points of the chains at time t, measured by the 
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Figure 3: (Color Online) The fidelity (F) between the the output and the input states 
averaged over all input states, as a function of kT and time (t), in a fixed magnetic 
field (B = 1), for an anti-ferromagnetic channel (J = —1). 

concurrence of the density matrix po',3- & fixed magnetic field, the concurrence of 
this density matrix is a function of temperature. Figures @ and © show this 
concurrence as a function of temperature and time for the ferromagnetic and anti- 
ferromagnetic chains respectively. All the previous comments apply also to this type 
of behavior. The striking difference is that in the anti-ferromagnetic chain, there are 
long intervals of time when no entanglement can be distributed in the chain regardless 
of the temperature, entanglement transfer is possible only in short periods of time. 
In fact comparison of the figures for the fidelity and concurrence shows that when 
the fidelity of the channel drops below the approximate value of 0.6, it no longer can 
transfer any entanglement. 

5 Summary 

We have studied the effect of thermal fluctuations on a recently proposed method for 
transportation of unknown states through quantum spin chains. We have developed 
a low temperature expansion which can be used to calculate this effect to a desired 
degree of accuracy at any given temperature. As an example we have calculated 
exactly the effect of thermal fluctuations on transportation of states on a short spin 
chain and have shown that the optimal time of extraction of transported states at the 
end of the chain is almost independent of temperature, the only effect of which is to 
lower slightly the fidelity of the output state with the input state. We have made a 
detailed comparison between the ferromagnetic and anti-ferromagnetic channels. 
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Figure 4: (Color Online) The concurrence C of the state of the endpoints of the chain 
as a function of time t and kT in a fixed magnetic field for a ferromagnetic chain. 
Here B = 1 and J = 1. 




Figure 5: (Color Online)The concurrence C of the state of the endpoints of the 
chain as a function of time t and kT, in a fixed magnetic field {B = 1), for a anti- 
ferromagnetic chain (J = —1). 
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7 Appendix : Spectrum of three and four Site Spin chain 

In this appendix we collect the eigenstates and eigenvalues of the hamiltonians He 
and H shown in equations ()31() and H32|) . 

The eigenstates and eigenenergies of the channel are as follows: 

|a2) = -(!+,-,-) -\/2|-, +,-) + !-,-,+), 
I«3) = ^{\+,-,-) + V2\-,+,-) + \-,-,+), 
l«4> = -^(1+,-,-)-!-,-,+)), (35) 



and 



with energies 



\ai+i) = af\ai) i = !,•••, 4, (36) 



El 
E2 



J 3B 

2 ~ ~Y' 
J B 



2 2 ' 

E> = -f , (37) 

and 

E,+4{J,B) = E,{J,-B) i = !,•••, 4. (38) 

The eigenstates of the total hamiltonian H (equation H32|)) are obtained by using the 
symmetries (|33|) . We use the notation |i) or \i,j) to indicate that the spins in the i— 
th position or the (i,j) positions are up and the rest are down: 



xi) = |-, -,-,-) 

X2) = ^(|1) + |2) + |3) + |4)) 

Xs) = ^(|1)-|2)-|3) + |4)) 

X4) = ^-^ (|l)-(^/2 + l)(|2)-|3))-|4)) 
2V2 + ^2 
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1X5) = 2_^ (|l) + (^-l)(|2)-|3))-|4)) 
1X6) = -^(|1,4)-|2,3)) 

1X7) = -^(|1,2) + |1,3) + |1,4) + |2,3) + |2,4) + |3,4)) 

1X8) = ^ (|1,2) - (1 + V2)(|l,3) - |2,4)) - |3,4)) 

2V(2 + \/2) 



1X9) = , ^ (-|1,2) + (1-V2)(|1,3)-|2,4)) + |3,4)) 



2^(2-^/2) 

Ixio) = ^-^(|l,2)-7?+|l,3)+e+|l,4)+e+|2,3)-r/+|2,4) + |3,4)) 

Ixii) = ^-^(|l,2)-7?_|l,3)+e+|l,4)+e-|2,3) + -7?_|2,4) + |3,4)) (39) 

where ^-t := 1 it \/3 and = 2 it \/?>. The other five states are obtained by the 
action of the flip operator a^^'^ on the first five states above, that is: 

IXi+ii) = Ixi) i = !,•••, 5. (40) 
The energies of the above states are: 

Ex = -^J-2B E2 = -jJ-B E3 = ^J-B 

E4 = ^{1 + 2V2)J-B E5 = ^{1-2V2)J - B Eq = 

E7 = -^J Es = ^{l + 2V2)J Eg = ^{l-2V2)J 

Eio = ^V+J E,x = -^r^_J. (41) 

and 

Ei+n{J,B) = Ei{J,-B) i = !,■■■, 5. (42) 
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